The Search for Cosmological Black Holes: 
A Surface Brightness Variability Test 
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ABSTRACT 

Recently it has been suggested that the majority of dark matter in the universe resides in the 
form of Jupiter mass black holes distributed cosmologically. This population makes itself apparent 
by microlensing high redshift quasars and introducing pronounced variability into their observed 
light curves. While several arguments dismissing this hypothesis have been presented, a conclusive 
observational test is, alas, sadly lacking. In this paper we investigate the effect of a cosmologically 
distributed population of microlensing masses on galaxies at low to intermediate redshift. The 
magnification of bright stars in these galaxies leads to small, but observable, fluctuations in their 
surface brightness. The variability time scale for Jupiter-mass lensing objects is of the order 
of a few months and this population may be detected through a future space-based monitoring 
campaign of a field containing z ^ 0.5 galaxies. The monitoring of galactic surface brightness 
will provide an effective test of the nature of dark matter on cosmological scales. 

Subject headings: Microlensing; Dark Matter: Black Holes: Galaxies 



1. Introduction 

The nature of dark matter is one of the ma- 
jor outstanding problems in current astrophysics, 
with the community being divided between two 
broad camps: those who propose particle dark 
matter and those who support massive compact 
objects. Both these ideas have received recent 
impetus due to the measurement of a neutrino 
mass (Fukuda et al. 1998) and the detection of 
MACHOs in the halo of our Galaxy (Alcock et 
al. 1993), although the small mass of the neutrino 
and the uncertainty in the location of the MACHO 
lenses suggests that the nature of majority of the 
dark matter still remains hidden from us. 
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One recent suggestion proposes that the vari- 
ability seen in high redshift quasars, rather than 
being intrinsic to the source, is due to the mi- 
crolensing action of Jupiter mass black holes dis- 
tributed cosmologically along the line-of-sight to 
the quasar (Hawkins 1993, 1996). To account 
for the characteristics of the variability, a sub- 
stantial population of black holes is required, 
with a density high enough to account for a sig- 
nificant fraction of the total mass in the uni- 
verse (Schneider 1993). While there is theo- 
retical support for such a model, with models 
of the quark-hadron transition producing black 
holes with a mass of ~ 10^'^ Mq [e.g. Hawkins 
& Taylor (1997)], observations have failed to re- 
veal their presence in studies of the equivalent 
widths of distant quasars (Canizares 1982; Dal- 
canton et al. 1994), although this conclusion has 
been contested (Hawkins 1996). Other objections 
to this microlensing hypothesis have been raised 
(Baganoff & Malkan 1995; Alexander 1995), with 
corresponding rebuttals (Hawkins & Taylor 1997). 
More recently, arguments have become quite vitri- 
olic with the suggestion that this idea of cosmolog- 
ically distributed black holes is being ignored by 
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the astronomical community cvcm though the ap- 
parent evidence is overwhelming (Hawkins 1998). 

The microlensing of cosmological supernovae 
provides a potential probe of any population of 
black holes (Schneider & Wagoner 1987; Ranch 
1991,a; Metcalf & Silk 1999), although their rarity 
suggests that a long observational period may be 
required until an unambiguous detection is made. 
In this paper we provide a complementary probe of 
a cosmologically distributed black hole population; 
instead of supernovae as the source population, we 
focus here on microlensing on the much more nu- 
merous population of normal stars in galaxies out 
to z ~ 0.5. In Section 2 we outline the analytic 
and numerical approach to this problem, while in 
Section 3 the results of a suite of simulations are 
presented. Section 4 is concerned with the time 
scale of the microlensing variability, and the ob- 
servational feasibility of this study is discussed in 
Section 5. Section 7 presents the overall conclu- 
sions of this study. 

2. Method 

2.1. Background 

Microlensing surveys of the Galactic halo focus 
upon identifying brightness fluctuations of a sin- 
gle, isolated source by an individual star or binary 
system along the line-of-sight (Paczynski 1986). It 
is not necessary, however, to identify single stars in 
a search for microlensing. Proposing to search for 
MACHOs in the halo of M31, Crotts (1992) noted 
that while the light registered in a single pixel is 
a combination of the entire stellar population, the 
majority of it arises in a relatively small number of 
extremely luminous giants. In a stochastic manner 
a massive object in the halo of M31 will traverse 
the line-of-sight to a giant and induce a strong 
microlensing magnification, increasing the surface 
brightness of the entire pixel. A monitoring pro- 
gram of this 'pixel lensing' has recently come to 
fruition with the identification of several poten- 
tial microlensing events (Crotts & Tomaney 1996; 
Tomancy & Crotts 1996; Alcock et al. 1999). 

Given the relative lens geometry, however, 
the optical depth for such lensing is very small 
(10~^ — > 10"'''), and events arc extremely rare. 
Similarly, if all dark matter is in the form 
of MACHO-like objects distributed uniformly 
through-out the universe, the additional optical 



depth these would add in our view to M31 would, 
due to the small distance and hence low density 
of material, be negligible. Focusing instead out to 
intermediate redshift, however, the optical depth 
becomes more substantial (this is demonstrated 
below) and a number of microlenses will lie along 
our line-of-sight to a distant galaxy. In such a 
regime a fraction of the population of stars in a 
pixel will be subject to strong magnification. As 
the intrinsic motions of the microlensing objects 
change their positions with respect to the line- 
of-sight, the 'state' of the microlensing (i.e. how 
much each star is magnified) will also change, re- 
sulting in a fluctuation in the apparent luminosity 
of a patch of stars. The details of these fluctu- 
ations depends on several factors, including the 
optical depth to microlensing and the luminosity 
function of stars in the source galaxy; these are 
explored in more detail below. 

2.2. Cosmological Considerations 

Standard cosmological models assume that all 

material in the universe is smoothly and homoge- 
neously distributed. In the following argument, 
however, all the matter in the universe is con- 
densed in to compact objects. To begin with, it 
is assumed that these objects do not influence our 
line-of-sight view to a distant source. Dyer and 
Roeder (1972;1973;1974) considered the propaga- 
tion of light through such an "empty beam" uni- 
verse and found that the convergence, or focusing, 
of the beam is removed due to the lack of material 
and that objects at a particular redshift appear 
fainter than those observed in a universe where 
the matter is smoothly distributed, even the uni- 
versal curvature is the same in both cases ^. For 
the purposes of this paper, we seek an intermedi- 
ate position, with matter in the form of compact 
objects distributed randomly in the universe, with 
no constraint on their location with respect to our 
line-of-sight to a distant source. With this, the mi- 
crolensing masses can influence the beam and, in 
the mean, the appearance of distant objects must 
be identical to that of a universe where all the 
matter is smoothly distributed (Weinberg 1976). 
With respect to the empty-beam universe model, 
the influence of the lensing masses is to magnify 

^An excellent description of cosmological distances is given 
by Schneider, Ehlers & Falco (1992) 
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the distant sources by a factor 

"ri(z) 



r{z) 



(1) 



where r{z) is the normalized angular diameter dis- 
tance to a source at redshift in a universe with 
a smooth distribution of matter, and ri{z) is the 
equivalent distance in a universe with the same 
global value of fl, but with all the matter in the 
form of compact objects that do not influence the 
line-of-sight view to distant objects. These angu- 
lar diameter distances are normalized such that 
D{z) = (c/HM^)- 

While the smooth matter case is analytically 
tractable, the solution of ri{z) typically requires 
the numerical solution of a differential equation 
(Dyer & Roedcr 1972; 1973; 1974). In the follow- 
ing, however, we will assume that Q = 1 and that 
all the matter is in the form of substellar com- 
pact objects. With these assumptions (Paczynski 
& Wambsganss 1989); 



ri{z) 3 + 2;2 + 2yrT^ + z(3 + VTTI) 



(2) 



r{z) 5(1 + z) 

While this tells us how bright objects will ap- 
pear in the mean, the distribution of microlcns- 
ing masses along any random line-of-sight will be 
different. This introduces a scatter in the bright- 
ness distribution and, as the relative positions of 
the microlensing masses and the source will change 
due to random motions of the lenses and sources, 
the magnification along any particular line-of-sight 
will vary with time. The first task is to determine 
the distribution of magnifications for a population 
of microlensing masses. 

2.3. Microlensing 

2.3.1. Microlensing Cross-Section 

To determine the degree of microlensing- 
induced fluctuations on the surface brightness dis- 
tribution of a distant galaxy, it is necessary to 
first derive the magnification (/i) probability dis- 
tribution for a source seen through a population 
of microlensing masses. This task is simplified 
by the fact that, irrespective of the mass spec- 
trum of the microlenses, in the regime where /i is 
large, p(/x) cx ^. In this case, the distribution 



^Additional caustic features do modify the magnification 
probability distribution in this region, although these fea- 



is dominated by bright pairs of images and can 
be treated analytically [e.g. Schneider & Weiss 
(1988)]. At lower magnifications, however, the 
form of the magnification distribution function is 
strongly dependent upon the optical depth to mi- 
crolensing, displaying complex secondary features 
due to higher order caustics (Rauch, Mao, Wambs- 
ganss & Paczynski 1992; Wambsganss 1992; Lewis 
& Irwin 1995); no analytic treatment fully char- 
acterizing these low magnification features has 
been presented [although recent semi-analytical 
attempts have proved promising e.g. Kofman et 
al. (1997) & Lee et al. (1997)]. 

For the case of a single, isolated microlens mass 
the probability distribution of magnifications can 
be derived analytically (Schneider, Ehlers & Falco 
1992) and is given by 



dfj, 



(3) 



Note that this probability is independent of the 
mass of the lensing object. It is assumed that, in 
the regime relevant to this work, the microlens- 
ing optical depth is small, and the masses do not 
influence one another to produce complex caus- 
tic patterns. With these assumptions, the global 
magnification probability distribution is simply re- 
lated to that of the isolated lensing mass. 

2.4. Normalization 

While Equation 3 presents the form of the mag- 
nification probability distribution, it still needs to 
be normalized before it can employed in any anal- 
ysis. This question was addressed by Linder, Wag- 
oner & Schneider (1988) who determined that in a 
universe with a global density f2, of which a frac- 
tion 1 — a is in the form of compact objects dis- 
tributed randomly with a constant comoving den- 
sity, then 



p{H, z) = 



-f2(l-a) 

r{z,Zs)r{zs) 
r{z) 



dz 



jl + z) 
^l + flz 



ri{zs) 



dm v{m)ps{fJ.,m) , (4) 



where i'{m) is the mass spectrum of the lensing 
masses. The last integration in this expression re- 
flects the fact that the maximum magnification a 



tures are negligible at the optical depths considered in this 
paper, a point we return to later. 
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source can undergo is related to its physical size, 
the mass of the lensing objects and the lensing ge- 
ometry (Chang & Refsdal 1979). Here, r{zi,Z2) 
is the normalized angular diameter distance be- 
tween two redshifts. For the gravitational lensing 
configurations presented in this paper, assuming a 
microlcnsing mass of lO"'' Mq and considering the 
analysis of Chang (1984), jjmax ~ 100 200/ VR, 
where R is the radius of the source star in So- 
lar radii. The most luminous B and A stars have 
R ~ lOi?0, for which Umax ~ 30 ^ 60, while 
the stars contributing the majority of the popu- 
lations luminosity possess R ^ 1 ^ IORq. This 
upper limit of the magnification cuts off the proba- 
bility distribution function, slightly enhancing the 
distribution below fimax [e.g. Figure 4 of Lewis 
and Irwin (1995)]. As significant magnification of 
stars can occur, we chose, therefore, not to re- 
strict the probability distribution function. With 
this, Equation 3 is independent of the microlens- 
ing mass function, i/{m) can be replaced with an 
arbitrary S function and Equation 4 can be used 
to derive a normalization, k, of the magnification 
probability distribution. 

While Equation 4 provides a normalization, the 
magnification probability distribution as described 
by Equation 3 diverges as /x 1; this is be- 
cause in determining Equation 3 it is assumed 
that the lensing mass lies within an infinite plane, 
with which any finite impact parameter results in 
fi > 1. An infinite impact parameter is required 
for fi = 1. As we are limited to a finite amount 
of sky we require that p{fi) ^ as /i ^ 1. We 
choose, therefore, to modify Equation 3 with a 
function to remove the divergence at /U = 1. The 
form employed is 

p(M) = (l-e^-(^-^))^ / , (5) 
(m^ - 1) = 

where both the values /Xu, and k are functions of 
redshift. This Equation is subject to two con- 
straints, the first being the usual normalization 
of the probability distribution function, and the 
second being flux conservation between full and 
empty beam universes (Equation 1): 

/oo 
Pin) dfx = 1 (6) 

/OO 
M Pin) dn = {nth) ■ (7) 



However, as Equation 4 normalizes the asymptotic 
form of the probability distribution, the single free 
parameter in Equation 5, namely ^w, can be fixed 
with one of these equations, while the other can 
be employed to ensure consistency. 

The procedure employed in this paper is to use 
Equation 6 to determine Hw- A series of these 
probability distributions, for source redshifts of 
z = 0.1, 0.2, 0.3, 0.4 & 0.5, with Q = 1 all in the 
form of compact objects, are presented in Figure 1, 
while Figure 2 presents k and over the redshift 
range z = 0.05 0.5. All of the magnification 
probability distributions possess a sharp peak near 
/U = 1, coupled with an asymptotic pin) oc 
tail at high magnifications. These normalized dis- 
tributions were then used to determine (^). This 
mean magnification is compared to the theoreti- 
cally expected value (Equations 1 & 2) in Figure 3; 
there is an excellent agreement between the two 
approaches. This is due to the low optical depth 
(and hence simple) regime in which the analysis is 
undertaken, as at higher stellar densities, such as 
those responsible for the microlensing of multiply 
imaged quasars, individual stars no longer behave 
like isolated lenses and their combined effects can 
lead to extremely complex magnification patterns 
on a high redshift source (Wambsganss, Paczynski 
& Katz 1990). 

Considering the theoretically expected magni- 
fication it is simple to determine that the opti- 
cal depth to microlensing out to the redshift of 
interest. This ranges from ~ 0.002 a.t z = 
0.1 to ~ 0.04 at z = 0.5. These numbers 
are several orders of magnitude greater than the 
K* ^ 10^^ 10^^ expected of gravitational mi- 
crolensing in the Galactic halo and towards M31 
[c.f. (Paczynski 1986)], and hence, in this oxtra- 
galactic scenario, microlensing effects should be 
more apparent. While these arguments demon- 
strate the efficacy of microlensing, deriving the 
detailed behavior requires knowledge of the un- 
derlying stellar luminosity function as the number 
density of potential sources will dictate the level 
of observed ffuctuations in the stellar sample as a 
whole, a topic which we turn to in the next section. 

How does Equation 5 compare to previous stud- 
ies of the magnification probability distribution of 
microlenses distributed in three dimensions? The 
numerical studies undertaken by Rauch (1991 ;a) 
focus on substantially higher optical depths than 
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those considered in this paper, ahhough one case, 
with = 0.01 [Figure 2 in Ranch (1991)], can 
be compared to distribution presented in Figure 1 
for z = 0.2 (k* = 0.008). An excellent corre- 
spondence between the two distributions can be 
seen over the entire magnification range. Kofman 
et al. (1997) and Lee et al. (1997) analytically 
tackled the form of the magnification probability 
distribution for both 2-D and 3-D distributions of 
microlensing masses. Also considering the numer- 
ical simulations of Rauch (1991; a) and Ranch et 
al. (1992), these papers investigate the nature 
of the higher order caustic features identified in 
the magnification probability distribution at large 
magnifications, features which are not accounted 
for in Equation 5. While substantially infliicncing 
the form of the 2-D distributions, these features 
were found to be very small for the corresponding 
3-D distributions. The analysis of Kofman et al. 
(1997) and Lee et al. (1997) was not considered at 
the optical depths presented in this paper, but for 
= 0.1 the caustic features induce a deviation 
of only ~ 15% at /i ~ 30, where the distribution 
has fallen to 10~^ that of its peak value. At the 
optical depths considered in this paper, such devi- 
ations will be even more inconsequential, and can 
be ignored. Lee et al. (1997) comment on diffi- 
culty in determining the form of the magnification 
probability distribution at low magnifications, es- 
pecially in the case where the microlensing masses 
are distributed in three dimensions. Comparing 
Figure 6 of Kofman et al. (1997) with Figure 13 
of Lee et al. (1997), it can be seen that at small 
magnifications (logioip,) < 0.1) that the proba- 
bility distributions for the 2-D and 3-D case are 
very similar. In this regime, Kofman et al. (1997) 
provide a semi-analytic formalism which approx- 
imates the form of the distribution. Comparing 
this to Equation 5 for the z = 0.5, = 0.04 case 
presented in this paper, a maximum difference of 
~ 7% at log\o{ij) ~ 0.005. Hence Equation 5 has 
an excellent correspondence with previous results. 

2.5. Stellar Luminosity Function 

With current instrumentation, the angular di- 
ameter of galaxies at redshifts of z ~ 0.5 ensures 
that they will extend over only a small number of 
pixels, and the light recorded in each pixel is the 
sum of the flux from a large number of stars. Each 
star in the pixel will be microlensed by some de- 



gree, and examining Figure 1, most stars will suffer 
a magnification ~ 1, although a fraction of the 
population will be magnified by a much more sub- 
stantial factor. On average, however, the overall 
flux in a pixel will be magnified by (/Utft,), although 
at any instant the brightness of a pixel will de- 
pend on the relative positions of the microlenses, 
as this dictates how much each star in the pixel is 
magnified. These will change due to the intrinsic 
velocity dispersion of the microlensing bodies, re- 
sulting in fiuctuations in the observed fiux within a 
pixel about this mean value. To be a useful probe 
of the intervening matter, however, these fluctua- 
tions must exceed some observable threshold. To 
properly determine this, the luminosity function 
of the source stellar population needs to be con- 
sidered. 

If the luminosity function is comprised of only 
relatively low luminosity stars then a large num- 
ber is required to account for total flux detected 
in a pixel. When such a population is microlensed 
many stars will be in a highly magnified state. 
As the microlensing state changes, and the mag- 
nification of each star also changes, the absolute 
number of stars in the highly magnified state re- 
mains essentially constant. Even if a single star is 
extremely magnified, a very rare occurrence, then 
the contribution this star makes to the flux reg- 
istered in a pixel is small compared to the over- 
whelming background of the numerous stars which 
are only moderately magnified. Conversely, if the 
luminosity function consists of only very luminous 
stars then only a small number are required to 
account for the luminosity in a pixel. It is very 
unlikely that any individual star is strongly mag- 
nified. When this rare event does occur, however, 
the contribution from the magnified star can dom- 
inate the flux received in a pixel, resulting in a sub- 
stantial fluctuation in the observed surface bright- 
ness. The real situation, however, is a combination 
of the above and we expect that a stellar popula- 
tion sampled in a pixel consists of a large pro- 
portion of faint stars which are insensitive to the 
effects of microlensing, resulting in a background 
flux, and a relatively small number of very lumi- 
nous stars which dominate the luminosity. 

The luminosity function model we choose is the 
stellar LF from Jahreil3 and Wielen (1997), derived 
from Hipparchos data within 20 pc of the Sun. 
The characteristics of this luminosity function are 
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presented in Figure 4; representing a population 
with a total luminosity of 10^ L©, the solid line 
shows the number of stars possessing a luminosity 
greater than some value, while the dot-dashed line 
is their contribution to the total luminosity of the 
pixel. It is obvious from this picture that while 
stars with L > 10 L© represent less than 1% of 
the population by number, they account for 70% 
of the total luminosity. 

3. Magnification Distributions 

Given a population of microlensing masses, the 
magnification of a single star can be chosen by se- 
lecting from the magnification probability distri- 
bution given by Equation 5. For a sample of stars, 
magnifications can be drawn from this distribution 
multiple times to determine the total magnifica- 
tion of the population. While this is straight for- 
ward, it becomes computationally expensive when 
the population consists of a very large number 
stars. The procedure employed in this study is 
to determine the magnification statistics for vari- 
ous numbers of source stars, combining the residts 
to reproduce the effect of microlensing the entire 
source stellar population. Analytically this can be 
calculated by convolving the magnification prob- 
ability distribution for a single star (Equation 5) 
with itself Usamp times, where Ugamp is the num- 
ber of stars being microlensed. Again, while this 
is simple for a small population of source stars, 
the calculation becomes numerically unwieldy for 
a large number of stars. For the purpose of this pa- 
per, these statistics were generated a Monte Carlo 
approach, drawing Ugamp observations from Equa- 
tion 5 and combining them to determine the mean 
magnification of the sample. Repeating this proce- 
dure 10^ times for each Usamp it is then possible to 
determine the distribution of mean magnifications. 
The results of this sampling for several source red- 
shifts between z = 0.05 and z = 0.5 are presented 
in Figure 5, with sample sizes ranging from 1 to 
10^ source stars. At all redshifts the same trend 
is seen with the distribution of the means evolv- 
ing from the single realization for (Equation 5) to 
a more Gaussian form as the sample size is in- 
creased; this is a consequence of the central limit 
theorem. These distributions can then be com- 
bined to determine the distributions of means for 
any sample size. 



Armed with this information it is now possi- 
ble to tackle the question of the overall bright- 
ness fluctuations expected for a population of stars 
drawn from a luminosity function, and hence the 
expected surface brightness fluctuation within a 
pixel. For this, the above magniflcation distri- 
butions are convolved with the stellar luminosity 
function described in Section 2.5, assuming a total 
luminosity, and hence number of stars in the pop- 
ulation. Again, a Monte Carlo approach is taken 
and the luminosity function is binned logarithmi- 
cally. While the bins at lower luminosity contain 
many stars they contribute a negligible fraction of 
the total luminosity and are uniformly magnifled 
by (/i) . At the high luminosity end, where the lu- 
minosity contributions of the bins become a more 
appreciable fraction of the total luminosity, magni- 
flcation probability distribution for each bin is cal- 
culated by combining the appropriate population 
magnification probability distributions presented 
in Figure 5. The distribution for each bin is then 
combined, weighted by the luminosity fraction of 
the bin, to determine the magnification probabil- 
ity distribution for the entire population of stars. 
Two constraints, that the mean value of the final 
distribution must be (n) and that the integrated 
area under the resultant probability distribution 
must be unity, are checked to ensure consistency. 
Figure 6 presents the results of this procedure for 
a range of source redshifts between z = 0.05 and 
z = 0.5, displaying the percentage variability for 
stellar populations of 10^ L© (dark grey), 10^ Lq 
(medium grey) and 10^ L© (light grey) total lumi- 
nosity. Table 1 presents the integral properties of 
these distributions with each column representing 
the range within which 10, 30, 50, 70 & 90% of 
the magnifications of the population lie, relative 
to the average magnification. A number of inter- 
esting features are immediately apparent; at low 
redshift the relatively small microlensing optical 
depth induces negligible fiuctuations into the ob- 
served total brightness of any of the populations. 
This changes as the redshift, and hence the mi- 
crolensing optical depth, is increased, until by red- 
shift z = 0.5 microlensing will introduce fluctua- 
tions of order ^ 1% into a sample with a total 
luminosity of 10^ Lq. Increasing the luminosity of 
the stellar sample, however, reduces the width of 
the distribution, as expected from the central limit 
theorem, such that a population of 10'' at z = 0.5 
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will suffer fluctuations of ^ 0.5%. Wc turn to the 
observational aspects of such small variations in 
Section 5. 

4. Event Timescales 

While the magnification probability distribu- 
tion is independent of the form of the underly- 
ing mass spectrum of lensing objects, the time- 
scale of any fluctuation depends implicitly on 
the mass of the microlenses. At the low opti- 
cal depths considered here, any variability will be 
very simple in form, consisting of the superpo- 
sition of isolated peaks rather than the complex 
light curves expected for high optical depth mi- 
crolensing [e.g. Lewis, Miralda-Escude, Richard- 
son & Wambsganss (1993)]. As with the mi- 
crolensing of Magellanic Cloud sources, we can de- 
fine a variability timc-scalc to be: 



(1 + ^/) 



4GM DuD,i 



Do 



23.4 



(1 + ^0 

1^1000 



M 
Mo 



r{zi, Zs)r{zi) 
r{zs) 



h- 



where M is the microlensing mass, iiiooo is the 
relative transverse velocity in units of 1000 kms"""^, 
and r(zi, Z2) is the angular diameter distance. 

For source galaxies out to 2; ~ 0.5, and consider- 
ing a lO"'' Mq microlensing mass, this timc-scalc 
is presented in Figure 7. The time-scale of the 
microlensing events are substantially less than a 
year for all lensing configurations and hence can 
be observed in a single observing season. When 
examining Figure 7, it is important to remember 
that the distribution of microlenses arc volume- 
weighted (c.f. Equation 4), and hence are more 
likely to be located nearer to the source redshift. 
Similarly the typical timc-scalc for the microlens- 
ing events seen in the light curve of a population 
of stars will also be skewed to values less than the 
maximum seen in Figure 7. Detailed light curves 
are beyond the scope of this current paper and 
they will be the subject of a further study of this 
phenomenon. 

5. Observational Considerations 
5.1. Variability Detection 

While the previous theoretical analysis has pre- 
sented the framework for determining the effects of 



a cosmological population of microlensing masses, 
it is important to consider whether the amplitude 
of the expected variations are observable. Here we 
consider monitoring many distant galaxies simul- 
taneously with a CCD detector, with resolution 
elements that each cover a population of luminos- 
ity - 10'^ L0. I gnoring K-corrcctions, this corre- 
sponds to an apparent magnitude of m = 29.0 at 
z = 0.5 (assuming qo = 0.5 and h = 0.75). It is 
immediately clear that such observations will nor- 
mally be sky-limited, except for the Next Gener- 
ation Space Telescope (NGST) for which the sky 
background is expected to be very low, and the 
spatial resolution high enough to have very small 
pixels. 

To detect 1% variations at the 3cr level, 0.33% 
relative photometry is needed, that is, S/N = 333. 
Assuming an 8m diameter mirror, and a diffrac- 
tion limited system, a count rate of ~ l.SSej / pho- 
tons/sec is expected from a m = 29.0 object in 
the J-band with NGST, where Seff is the to- 
tal system efficiency. For exposure times longer 
yr{^] ^^^^ ^ f^'^ minutes, such observations will be 
photon noise dominated, since a background of 
< 10~^ photons/sec/resolution element ^ is ex- 
pected at l.l^TO. So the required S/N is achieved 
in n /Seff hrs of exposure. For a population of 
stars at a redshift oi z ^ 0.5 with a total lumi- 
nosity of 10 ''L0, the probability of the necessary 
1% variation is 4%, so several thousand such pop- 
ulations need to be monitored to detect a sample 
of microlensing events. (Although deeper observa- 
tions, focusing upon the fainter regions of galaxies, 
will reveal similar scale fluctuations in ~ 23% and 
~ 56% of the pixels for populations of, respec- 
tively, 10® L0 and 10^ Lq). Of course, in a typ- 
ical field, many resolution elements on the CCD 
camera will cover z ~ 0.5 galaxies, and so the 
monitoring of the populations may be achieved si- 
multaneously. If NGST is resolution limited, each 
resolution element will cover 0.001 □". For each 
resolution element to cover a 1O''L0 population 
at z ^ 0.5, requires the observation of a region of 
surface brightness 21.5mag/n". In regions fainter 
than 21.5mag/n", boxes larger than a resolution 
element can be summed up, to give a total lu- 
minosity 10^ Lq; the variability of these boxes can 
also be studied with the penalty of a small amount 



^see http;//www.ngst. stsci.edu/sky/sky.html 
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of extra sky noise. 

The HDF field shows that the area in a typ- 
ical 2' X 2' field covered by surface brightness 
> 21.5mag/n" regions of z — 0.5 galaxies ex- 
ceeds several tens of □", that is, tens of thou- 
sands of NGST resolution elements. The above 
cosmological microlensing-induccd variability will 
be observed over a background of intrinsic vari- 
ability events, due to supernovae, variable stars 
and galaxy self-lensing, as well as due to the in- 
evitable Poisson noise in the observed counts. Un- 
like current studies of 'surface brightness fluctua- 
tions' method for distance determinations out to 
galaxies at - lOOMpc (Tonry & Schneider 1988; 
Thomsen, Baum, Hammergren & Worthey 1997; 
Lauer, et al. 1998) which essentially requires sin- 
gle epoch observations, to identify microlensing 
induced variability a monitoring program is re- 
quired, allowing the identification of supernovae 
by their light-curves and peak luminosities. Su- 
pernovae are also rare, as are variable stars with a 
luminosity sufficient enough to influence the total 
luminosity of populations L ~ 10^1/0, if a sur- 
vey was to focus on more luminous source pix- 
els. As with supernovae, these would also reveal 
themselves via their characteristic light curves. If 
the depth of the survey is increased such that 
L lO^L© pixels are probed, the ubiquity of 
their variations would rule against variations in 
the source population. In the case of self-lensing of 
galaxy stars by MACHOs in the halo of the same 
galaxy provides a negligible contribution as the 
optical depth is many orders of magnitude smaller 
than that due to the cosmological microlenses con- 
sidered here. The main source of detected back- 
ground events is likely to be due to simple Poisson 
photon noise. However, in the situation outlined 
above, Scr variations due to Poisson noise are ap- 
proximately 40 times less common than the Scr mi- 
crolensing variations. Data at two epochs would 
allow an initial feasibility study. If the required 
variability rate is indeed observed, follow-up ob- 
servations at further epochs should be taken to 
monitor the light curves of the variability events, 
to distinguish microlensing events from unlensed 
supernova or nova events, and to reject variations 
due to Poisson noise. 



5.2. Extreme Events 

The previous Section has demonstrated that 
the typical fluctuations in the brightness of a pixel 
is detectable and lends itself to future space-based 
observatories. But are there any aspects of mi- 
crolensing by a cosmologically distributed popu- 
lation that would make them apparent with cur- 
rent technology? As noted earlier in this paper, 
while the vast majority of stars will be magni- 
fled by a value close to the theoretically expected 
mean (Equation 1), very rarely a star will un- 
dergo an extreme magnification. The probabil- 
ity that a particular star is magnified by an ex- 
treme value is found by integrating the magnifica- 
tion probability distribution (Equation 5), which 
is presented graphically in Figure 8. For instance, 
in a 10''' L0 population, the LF of Figure 4 gives 
~ 0.1 stars with L ~ 10^ Lq; if such a star were 
to be magnified by a factor of 20, it would alter 
the surface brightness of the population by 200%. 
At a redshift of z = 0.05, this is an uncommon 
phenomenon, with probability P(> = 10~^, 
however, the apparent magnitude of the popula- 
tion is much brighter: m = 27.5, and the pho- 
tometric accuracy needed to resolve these varia- 
tions is S/N ^ 3. While these are within the range 
of HST's capabilities, and many tens of fields of 
depth of m = 27.5 with S/N ^ 3 have been reim- 
aged, it is unlikely that a sufficient area has been 
covered to detect a sample of such rare events. So, 
while the tools are available, more substantial deep 
survey areas are required in a search for extreme 
events with HST. 

6. Other Cosmological Models 

While throughout this paper a standard flat 
cosmology with flo = I has been employed, re- 
cent studies of high redshift supernovae suggest 
the presence of a substantial cosmological con- 
stant (Garnavich, et al. 1998). To investigate 
the influence of such a dominant cosmological 
constant on the result presented in this paper 
the formalism of light propagation in generalized 
Priedmann universes was employed (Linder 1988). 
The full and empty beam distances in three uni- 
verses were calculated in three cosmological mod- 
els, (Qo = l,Ao = 0), (fio = 0.3, Ao = 0) and 
(fio = 0.3, Ao = 0.7), used to calculate the mean 
relative magnification between the full and empty 
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beam cases (Equation 1) and the niicrolcnsing op- 
tical depth. These are presented in Figure 9. 

Out of these three models, the case presented 
in this paper consistently displays higher optical 
depth over the entire redshift range; this is simply 
due to the fact that in this model more material 
lies along the line of sight to a distant source. In 
comparison, the optical depth in the open case 
is several times smaller, reaching a ~ 0.15 at 
z = 0.5. The addition of a dominant cosmolog- 
ical constant into such a universe does not greatly 
change the values of the optical depths. This is 
because, while the cosmological constant changes 
the global curvature of the Universe, it has no in- 
fluence on the Ricci focusing of the beam and lo- 
cally the Universe appears as it is solely matter- 
dominated (it is for this reason that cosmologi- 
cal supernovae need to be identified at redshifts 
greater than 0.5 in an attempt to determine the 
geometry of the Universe). 

7. Conclusions 

Recent studies have claimed that the fluctua- 
tions observed in the light curves of high redshift 
quasars are due to the microlensing effect of a 
population of cosmologically distributed Jupiter 
mass black holes. While arguments have been 
presented to refute this claim, no observational 
test of this hypothesis has been readily available. 
In this paper we have investigated the effect this 
putative population would have on the surface 
brightness distribution of galaxies out to interme- 
diate redshift. The results of this study conchi- 
sively demonstrate that if a significant population 
of cosmologically distributed compact objects are 
present, with enough mass to account for the uni- 
versal dark matter budget, then they will induce 
a 'twinkling' into the observed surface brightness 
distributions of galaxies at low to intermediate 
redshifts. The magnitude of the observed fluctu- 
ations are a function of the redshift of the source 
galaxy, ranging from negligible values in the local 
Universe, to ~ 1% at z = 0.5, with a time-scale of 
variability on the order of weeks to months. 

While the optical depth to microlensing at such 
low to intermediate redshifts is still quite small 
and the induced fluctuations are at a relatively 
low level, the sky density of galaxies out to these 
redshifts vastly exceeds that of supernova and 



quasars, the current focus for the search for cosmo- 
logically distributed microlensing masses. With 
sufliciently deep imaging with source pixel lumi- 
nosities of greater than IO^Lq, microlensing in- 
duced surface brightness variability can be de- 
tected over galaxies at intermediate redshifts, and 
the ubiquitous nature of this variability means 
it will dominate over any contaminating effects, 
such as supernovae, variable stars and self-lensing. 
The scale of the induced variability make obser- 
vations conducive over a single observing season, 
and with the advent of the Next Generation Space 
Telescope, such observations will soon be possible. 
Hence, a monitoring program to search for fluctu- 
ations of the surface brightness of the plethora of 
galaxies at z ^ 0.5 offers a effective test of the exis- 
tence of cosmologically distributed compact mat- 
ter. 
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Table 1 









10% 


30% 


50% 


70% 


90% 






lO^Lo 


±0.010 


±0.027 


±0.040 


±0.052 


±0.065 




=0.05 


IQf'Lo 


±0.004 


±0.013 


±0.021 


±0.029 


±0.041 








±0.002 


±0.005 


±0.010 


±0.014 


±0.021 








±0.021 


±0.063 


±0.101 


±0.138 


±0.186 


Z: 


=0.10 


IO^Lq 


±0.010 


±0.031 


±0.051 


±0.075 


±0.108 








±0.003 


±0.013 


±0.022 


±0.033 


±0.052 








±0.047 


±0.139 


±0.235 


±0.335 


±0.480 


Z: 


=0.20 


lO^Lo 


±0.022 


±0.070 


±0.122 


±0.177 


±0.266 






IO^Lq 


±0.009 


±0.027 


±0.046 


±0.071 


±0.116 






IO^Lq 


±0.072 


±0.215 


±0.365 


±0.530 


±0.786 


■L- 


=0.30 


lO^L© 


±0.036 


±0.104 


±0.187 


±0.277 


±0.427 






lO^Lo 


±0.010 


±0.040 


±0.070 


±0.108 


±0.175 






IO^Lq 


±0.092 


±0.285 


±0.490 


±0.720 


±1.094 


Z: 


=0.40 


lO^Lo 


±0.047 


±0.143 


±0.252 


±0.373 


±0.578 






IO^Lq 


±0.016 


±0.052 


±0.088 


±0.137 


±0.233 






IO^Lq 


±0.112 


±0.352 


±0.608 


±0.899 


±1.377 


Z: 


=0.50 


lO^L© 


±0.049 


±0.168 


±0.305 


±0.459 


±0.715 






lO^L© 


±0.005 


±0.057 


±0.108 


±0.159 


±0.279 



Note. — The integrated properties of the probabiUty distributions 
given in Figure 6. The cohimns presents the range over which 10, 30, 
50, 70 & 90% of the magnifications are found, relative to the mean 
value (i.e for a lO^L© population at z=0.40, 50% of all magnifications 
lie within ±0.252% of the mean value). 
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Fig. 3. — A comparison between the theoretical value of the image magnification, {nth), compared to that 
derived from the analytic approach presented in Section 2. 
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Fig. 4. — The properties of the assumed luminosity function, assuming a population with a total luminosity 
of IO^'Lq. The solid line represents the number of stars in the population with a luminosity greater that 
L*, while the dot-dash line is the summed luminosity of this subset of stars. While stars with L > IOLq 
represent only a small fraction of the over all population by number they are responsible for the majority of 
the total luminosity. 
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</x>=1.001 </2>=1.005 </2>=1.017 




</2>=1.034 </x>=1.056 <//> = !. 082 

Fig. 5. — The distribution of the mean magnification of a sample of stars for redshifts between z = 0.05 and 
z = 0.5. The dotted vertical line represents the mean value of the relevant magnification distribution, while 
the lower boimd on each plot is /i = 1. Each line consists of a stellar sample of 1, 10, 100, 1000 and 10000 
stars, depending on its shading. The scale factor, f, is introduced such that all the curves can be placed on 
the same y-scale. As these are probability distributions, the integrated area under the curve (accounting for 
the f-scaling) is unity. 
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Percentage Variability 

Fig. 6. — The distribution in total luminosity of stellar samples from 10^ Lq (broadest - dark grey) to 10^ Lq 
(narrowest - light grey), defined by the luminosity function in the text. At very low redshifts, the low 
microlensing optical depth induces only negligible variability, while by a redshift of 2: = 0.5, fluctuations of 
order a percent are induced. The integrated properties of these distributions are given in Table 1. 
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Microlens Redshift 

Fig. 7. — The variability time scale for a 10^'^ M© microlensing mass with a relative transverse velocity of 
1000kms-\ for sources at 2 = 0.1, 0.2, 0.3, 0.4 & 0.5 . 
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Fig. 8. — The cumulative magnification probabilities for sources at 2; = 0.05, 0.1, 0.2, 0.3, 0.4 & 0.5 . 
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Fig. 9. — The optical depth out to ^; = 0.5 in several cosmologies. At the low redshifts considered, the 
introduction of a cosmological constant has only a small influence on the optical depths. 
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